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ABSTRACT EGY 


A "split-field" infrared optical system has been designed and constructed to 
provide simultaneous image pairs in a single frame of an infrared (IR) imager, 
differing only in the direction of linear polarization. The optical train is afocal, 
allowing its use with a variety of infrared imaging devices. The system can operate 
in both long-wave IR ( 8-12 um) and mid-wave IR (3-5 um) with interchangeable 
polarizing splitter plates. Previous work at the Naval Academic Center for Infrared 
Technology (NACIT) at the Naval Postgraduate School has demonstrated that 
significant image improvement in infrared image contrast can be obtained by use 
of polarization filtering, especially for targets at sea through suppression of 
polarized sea background emission in the long wave, or of polarized reflection in 
the midwave. That work utilized digital subtraction of sequential image pairs with 
orthogonal polarizations, but suffered from inability to obtain simultaneous 
images and also from problems due to reflections from external polarization filters. 
Both of these problems are eliminated with the new split-field technique. Prelim- 
inary tests of this system with an AGA-780 imager were carried out in a field 
experiment using an Arleigh Burke DDG class ship as target at varying distances 
at sea. For comparison, images were also taken in successive pairs in time with the 
AGA-780 viewing the scene directly ( without the split-field adjunct) with 
interchangeable internal polarizing filters in a rotatable filter wheel. Subjective 
and numerical analysis of the data from the field experiment demonstrated good 


image quality and contrast improvement, and potential for future utility. 
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I. INTRODUCTION 


Previous work done at the Naval Postgraduate School has 
shown that infrared image contrast for sea targets can be 
improved by polarization filtering (Refs: 20,21,22). 
Unwanted background radiance which clutters the image can 
consist of emission and reflection components. At grazing 
angles an appreciable part of the water surface emission is 
linearly polarized with the electric vector vertical, while 
direct or indirect solar radiation reflected from the 
surface shows horizontal polarization. Other elements of the 
scene, such as the sky and man-made objects, often emit or 
reflect radiation that is unpolarized or polarized in a 
different direction from that of the sea surface. As a 
result, viewing the scene through a polarizing filter, or 
subtracting two otherwise identical images except for 
orthogonal polarization directions, can improve the target 
contrast against background. 

Polarization filters for the infrared are only recently 
avallable. These consist of "wire-grid" structures 
deposited by photolithography and etching on transparent 
substrates. Such substrates also need to be antireflection 
coated to reduce unpolarized reflection by the dielectric 


substrate. 


Special problems arise when dealing with IR as opposed 
to visible radiation. Because all materials emit radiation 
at temperatures above 0 degrees K, sensor background 
radiation often contaminates the final image. Very low 
temperatures must be used at all points from which radiation 
can enter the detector in order to reduce the amount of 
unwanted radiation entering below acceptable levels. Liquid 
Nitrogen and thermo-electrically cooled bodies are commonly 
utilized. Also, the 'narcissus effect’ can occur in the 
final image when the cold detector 'sees' a reflection of 
itself from a flat or back-focussing optical surface (Ref. 
23). Various methods to eliminate these problems are 
implemented in the techniques presented. These are key 
driving factors in the utilization of internal filtering and 
in the design of the Naval Postgraduate School Split-Field 
infrared polarized imaging system. 

Important to the study of image improvement is the 
quantitative characterization of image intensity 
differences. In previous experiments done at the Naval 
Postgraduate School in which the orientation of the 
polarizer had to be changed manually between recording 
orthogonally polarized images, there were changes in the 


image in the time interval during this process 


(Refs.1,21,22). Consequently it was impossible to have two 
identical images for comparison. To surmount this problem a 
Split-Image optical system was designed and built at NPS. In 
this system two orthogonally polarized images are recorded 
Simultaneously. The images can then be digitized and 
mathematically subtracted to demonstrate the contrast 
improvement. Specific design, implementation, performance 
and image analysis will be explained and presented. 

A description is also given of adaptation of the AGA 
imager with polarizing filters installed in the internal 
filter wheel. This system is useful in recording sequential 
images with orthogonal polarizations. 

This thesis explores the phenomenon of infrared 
polarized radiation, background and target polarization 
characteristics, and experimental techniques for filtering 


and improving image quality. 
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II. BASIC INFRA-RED THEORY 

To fully understand the phenomena associated with 
imaging of polarized infrared radiation it is necessary to 
understand the basic physics used to characterize infrared 
phenomena. This chapter deals with the laws of thermal 
radiation, radiation measurement metrics, and discussion of 
thermal imaging fundamentals. These short presentations 
attempt to give some insight into the selection of the 
infrared band and how it is used. 
7 INFRA-RED TERMINOLOGY AND UNITS 

In the study of radiation, the radiant power exchanged 
between the target, medium, and detector is of interest. A 
set of parameters is defined to describe the radiation at 
each point of interest. Table 2.1 shows some of the 
parameters and their units of measurement (Ref. 19). The 
last column in Table 2.1 denotes the units in the metric 
system. The entire content of this system of units will not 
be addressed, but some particular units need to be 
emphasized. Temperature is measured in Kelvins (K). The 


steradian (Sr) unit of solid angle is derived from defining 
the total solid angle around the center of a sphere as 4nz 


steradians. Note that these quantities are commonly used in 
inconsistent units. This is done to match imaging system 


dimensions more closely for system performance calculations. 









Symbol Name Description Units 
(metric) 
M Radiant Radiant Flux leaving Watts/cm* 
an element of a 
surface per unit 
(Exitance) area. 
E Irradiance | Radiant power per Watts/cm* 
unit area incident on 
iD Radiant Watts/Sr 
Intensity 


a surface. 
Beeoeee Watts/cm’sr 


Table 2.1 Radiometric Units 
All other physical parameters are measured in conventional 






















Radiant power leaving 
a point source per 
unit solid angle. 

















Radiant power leaving 
Or arriving ata 
surface at a point in 
a given direction per 
solid angle per area 
normal to that 
direction. 













units. Along the imaging path, M characterizes the power 


leaving the target surface per unit area. Some of this 


radiation is emitted by the target and some is reflected due 


to Irradiance from another source. The distant transport of 


the radiant energy from the target is described by the 
Radiant Intensity, I, which includes the geometrical 
spreading with range. JI is the most useful parameter in 


calculating the maximum range of an imaging detector. When 


the radiation reaches the detector, the Irradiance at the 
receiver of the target radiation 1S converted into a signal 
Lot TecoreinGg. 
B. LAWS OF THERMAL RADIATION 

In order to detect a target, it 1S necessary to 
characterize the radiation coming from the surface. Most 
infra-red detectors only operate efficiently in a small band 
of wavelengths, so the thermal radiation from a target must 
be known in that range. A set of well known radiation laws 
exists that describe the spectral emittance of thermal 
sources. 

as Black Body 

A black body is known simply as an ideal emitter and 


absorber. By Kirchoff's Law, for a black-body the spectral 
absorptivity(a) = emissivity(se) = 1, so the reflectance (p) 
and transmissivity (tT) must be = 0. For a black body this 


is true at all wavelengths. The spectral radiant emittance 
curve for a black body is shown by Fig. 1, and is described 


by Planck's Law. 


Spectral Radiant Exitance (w em? um”) 





0 5 10 15 20 25 30 
Wavelength (pm) 


Figure 1. Black Body Radiance 
2. Planck's Law 


Planck's daw 1s the distribution shown in Fig. a. im 


1S given by 


J Ee hye 
MAT) Saas (W/ cm*ym) (1) 


(e T_1) 


where h is Planck's constant, and k is Boltzman's constant. 
The spectral radiance of the black body can be reproduced by 


the Planck distribution. For non-black body emitters, the 


8 


radiant exitance is reduced by the emissivity(e). The 


radiant exitance is then 


MEM ack-body ( 2 


with e<xl for all non-black bodies. Examples of emissivities 














Table 2.2 Emissivities of common materials 


of practical materials are shown in Table 2.2 (Ref. 19) 
If ¢,(T) = constant independent of wavelength, then the body 
is called a grey body. For most surfaces this is not the 
case, and the spectral emissivity varies greatly with 
wavelength. Many materials behave as grey bodies over a 
limited range of wavelength. 

3. Wien Displacement Law 

In Fig. 1, a maximum radiant exitance occurs at a 


wavelength Amax- Amma, Can be found by differentiating the 


EPrenck GuscePOuUcion LUNCETensana equauing to zero.” Tne 


result is the Wien Displacement Law. 


A, T=2898.78pmK (3) 


Using Wien's Displacement Law, the wavelength with the 
highest spectral radiance is known with the surface 
temperature. 

4. Stefan-Boltzman Law 

If the radiant exitance is integrated over all 
wavelengths for a given temperature, the total power emitted 
per unit area can be found. The equation is found to be, 


frouyaubillack body, 
M(T)=oT' (4) 


o is the Stefan-Boltzman constant, o=5.67 x 10° W/cm*K’. 


For other bodies, application of the spectral emissivity 
must be done over all wavelengths. 
C. ELECTROMAGNETIC BANDS AND ATMOSPHERIC WINDOWS 

To understand why imaging is performed in the infra-red 
band, knowledge of the bands and the atmospheric attenuation 
"windows' is necessary. All bodies radiate thermally at all 
wavelengths, so it would appear logical to detect the 
radiation at the peak wavelength defined by Wien's 


Displacement Law. This is not done in every case, due 


LO 


partially to atmospheric effects on the emitted radiation. 
A brief discussion of how the target and atmosphere interact 
follows. 
1. Electromagnetic Radiation Bands 
The infrared region of the electromagnetic spectrum lies 
between the radio (longer wavelength) and the visible 
(shorter wavelength)regions. Fig. 2 shows this region of 


the spectrum (Ref. 19). Its bandwidth is approximately 


-lum to 1000um. The longest wavelength of interest for 


conventional IR imaging is approximately 15um since most 


targets and modern quantum detectors operate below this 
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Figure 2. Partial Electromagnetic Spectrum 
wavelength. termed the far infrared limit. The region 


nearest the visible is the near infrared. The "millimeter 


wave” band occurs around 1000um. 


dale 


2. Atmospheric 'Windows' 

The processes that effect the radiation on its path 
through the atmosphere are extremely complex and will only 
be considered briefly. Attenuation affects the amount of 
radiant power incident on the optical system, which in turn 
affects the overall performance of the system. The 
attenuation mechanisms of primary concern are absorption and 
scattering of the radiation. LIMLtmng Onpeavoiding the 
amount of absorption and scattering in the target path is 
essential for any thermal imaging system because of the 
severity of the effects. Suspended aerosol particles are 
the primary contributors to scattering of the radiation in 
the IR. Aerosols include dust, fog, salt, ash, and many 
other airborne particles. Molecules such as H,0 and CO, are 
the predominant absorbing species in the atmosphere (Ref. 
19) and define the windows in which IR radiation can pass 
through the atmosphere with any effectiveness. This is due 
to the structure of the molecules and the modes in which 
they can be excited and thus absorb energy. Figure 3 shows 
a plot of atmospheric transmittance versus wavelength over 
the near to far infra-red. This plot was generated by an 
atmospheric propagation modelling code, LOWTRAN 7. LOWTRAN 
7 has resident atmospheric models for aerosol and humidity 


content in addition many other atmospheric characteristics. 


| 7 


The model used for the plot is for a Mid-Latitude Summer 
atmosphere, Skm path at sea-level. The ‘'Rural' designation 
is for the aerosol model describing an air mass over land 
with primary components being dust, sand, and hygroscopic 


particles [19]. 
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Figure 3. Atmospheric Windows 


Transmission with no attenuation from the figure has a 
transmittance value of 1. The windows are depicted by the 
regions of maximum transmittance. It is clear that any 
imaging would need to take place for radiation found in 


these regions. Commonly used wavelengths for imaging 
military targets are around 4.0Oum (short wave) for hot 
sources such as exhaust plumes and power plant components, 
and 8-12um (long wave) for surfaces nearer to ambient 


temperatures. The relation of the temperature to desired 


wavelength band follows the Planck distribution for maximum 


iS 


radiated power. This correlation between radiated power 
wavelengths and propagation windows must be considered when 


using IR radiation for imaging is desired. 
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III. POLARIZATION PHYSICS 


Polarization occurs in the reflected and emitted IR 
radiation of targets and backgrounds. For purposes of 
enhancing target contrast (as exemplified by ship targets at 
sea), it is necessary to determine the degree of 
polarization present and its orientation. Fundamental 
physical laws and polarization equations describe the 
important parameters of polarization of the target, and sea 
and sky background. Quantification of these parameters 
leads to expected increases in contrast and is necessary in 
designing and implementing polarizing imaging systems. 

A. STOKES PARAMETERS 

Sir George Gabriel Stokes (1852) discovered that the 
polarization state of any light source could be completely 
described by four observable quantities which have become 
known as the Stokes Parameters. This is done by taking the 
time average of the E field vector amplitude described by 
the polarization ellipse, thus yielding an observable 


intensity (Ref.17). The resultant four parameters are 


eae Es 
& P 
M=E2-E? (3) 
C=2E E cosé 


S=2E ES ind 


H. and B are Che normal and parallel comeonents Of hem wae 


i 


electric field, with respect to the plane of incidence. 

For a horizontasl surface the vertical component (p) lies in 
the plane of incidence, also referred to as parallel 
orientation. JI denotes the total intensity of the light, mM, 
the intensity of linear 's' or 'p' polarized light, C, the 
linear +45 degree polarized light, and S, the circularly 
polarized intensity. 6 is the phase difference between the 
Ss and p components of the E field. The degree of 
polarization is the ratio of the intensities of the 


polarized to the total intensity [3]. 


e+ 2+g2) 1/2 
7 (6) 
I 


The polarization of the radiation to be filtered comes from 
effects due to emission and reflection. Specific laws have 
been developed to describe these phenomena and their 
contribution to the Stokes parameters will be shown. 
EB. EMISSION POLARIZATION 

The target and background radiance produced by thermal 
emission from surfaces is inherently polarized [9]. Of 
specific interest for imaging targets at sea is the emission 
polarization of the sea and sky background since the primary 
contrast improvement will be obtained by filtering of the 


background. 
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ae Sea Emission Polarization 
To find the emission polarization of the sea, the 
surface can be treated as an interface with coordinates 


shown in Figuse 4a(Refiee 21). 






Observer 


ware rere 


x 
Figure 4. Coordinate System of Sea Surface 


Y is the azimuthal direction of the sun. S and D are the 


directions of the sun and observer, respectively. N is the 
wave facet normal, @and@ the zenith and azimuth angles for 


each vector. The light coming from the upper side of the 
interface (+z) consists of a reflected and an emitted 


component, such that the total radiance is given by 


ile 


N(8) =N, (8) +N, (8) (7) 


The light emitted from the surface can be considered as 
radiation coming from the sea water side of the interface 
undergoing refraction at the surface. Behavior at the 
interface can be derived from the Fresnel equations and is 
discussed in this fashion by Sandus (Ref. 9) and Sidran 
(Ref. 11) following Reference 16. The emissivity of the 


surface is given by Reference 16. 


z, is dependent on the orientation of the polarization 


(parallel or perpendicular) where: 


eG 


stcos’O 
Zi= 


/2cos@[ (sta-sin’®) 1/7] 


st (a’*tb’*) cos’0 


a 
/2cos®@[s(a’7+b?)+(a-sin’®@) (a?-b?)+(2ab?) J” (9) 


a=Re(n7’) 


b=Im(n7) 
s=y (a-sin’@)*t+b° 


n is the complex index of refraction of the medium, ©, the 


angle of emission. These equations hold true for any 
specular surface. Dependence on wavelength comes from the 
index of refraction (Ref. 11). This analysis indicates that 
both the emitted and reflected components of surface 
radiance will show polarization, predominately in orthogonal 
directions. Note that consideration is given only to the 
linearly polarized light. For most materials, there is only 
a small range of angles at which elliptical polarization 
occurs. Even in these cases the phase difference between 
orthogonal components is small enough (due to small 
dielectric constants) that an approximation of linearly 


polarized light at all angles is valid (Ref. 3). This 


ne) 


theory has been applied to the sea surface and proven by 
Sidran (Ref. 11), Baesner and McCoyd (Ref. 8), and Gregoris, 
Yu, Goeper, and. Milne ma Reeweo).. 

In dealing with the sea surface, wave structure alters 
the coordinates of interest in Fig. 4. The wave structure 
1s considered to compose a distribution of emitting facets 
with unit direction toward the observer defined by the 


vector S (see Fig. 4), where 


S-s,its,jts,k (10) 


The orientation of the normals to the facets with respect to 


the normal to the plane defined by g and the z axis is 
defined by the angle a. The effective emissivities of each 
facet referred to the vertical and horizontal directions, 
are given by 


¢ CeO, sin’a 
Te 


€ : oS in Ot et cos’a 


The prime symbol denotes the emissivity in the observer 


direction relative to the vertical to the plane containing 
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the z axis and the observer direction g (Ref. 8). The 
distribution of the slopes of the facets has been 
experimentally determined by Cox and Munk (Ref. 15). The 
Gram-Charlier equation with empirically derived constants 
adequately describes the distribution of wave slopes (Ref. 


15). The resulting wave slope distributions are considered 
valid for wind speed <14m/sec and wave slope angle (f)<23 


degrees. Plots of the Cox and Munk distribution are shown 
in Fig. 5 from (Ref. 11). Because shadowing occurs on the 
back surfaces of the waves in the direction of the observer, 
a correction should be applied to the equations as used by 
Basener and McCoyd (Ref. 8) for percent polarization 
calculation. Saunders (Ref. 10) gives a correction to the 


Basener and McCoyd results. Figure 6 shows the polarization 


of emission for 1.0um wavelength and 14m/s wind speed after 





Sidran (Ref. 11). The percent of polarization, Q is 
defined by 
= ee, 
2,08 


iL 


The results in Fig. 6 hold true for all wavelengths < 10mm 


excluding the absorption bands (Ref. 11). Note that the 


Zl 


largest polarization effects are seen at grazing incidence. 

2. Sky Emission Polarization 

The IR radiation from the sky is unpolarized (Ref. 2). 
Use of polarized filtering components will not contribute to 
the contrast improvement measurements when the sky is viewed 
directly. The sky emission in the image will be subject to 
attenuation by the polarizing filter to the degree of the 
unpolarized transmissivity o£ the filter. Hag. 7 (Refi. ag) 
shows the sky radiance as a function of zenith angle. 
Specific radiance characteristics can be seen in Fig. 7. 
When viewed toward the zenith, the sky emits less due to the 
fact that the atmosphere consists of water vapor and 
aerosols (absorbers) which also act as good emitters, but is 


less dense. 
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Figure 7. Sky Radiance as a 
function of zenith angle from 
Saunders [10]. 
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Sky radiance increases with zenith angle toward the horizon 
to that of a black-body with atmospheric temperature. This 
is due to the maximum thickness of the atmosphere 
encountered at that angle and even the small emissivities 
of the absorbers combine (Ref. 10). 

Ze Target Emission Polarization 

Although primary contrast improvement in the at-sea 
scenario will be gained by filtering of sea emission 
polarization, there is polarization due to emission from the 
target itself that must be considered. Painted surfaces 
such as are found on all types of vehicles and ships display 
emission polarization (Ref. 3,4). Figure 8 shows 21 
measurements taken on various paint samples ranging from 


very smooth torough sand-paint mixtures (Ref. 4). 
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Figure 8. Ref [4] 
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The line for an ideal surface in Fig.8 denotes the amount of 
polarization that would be seen from a perfectly specular 
reflector. The rougher samples show less polarization than 
the smoother. The amount of polarization depends on the 
index of refraction of the paint and any degradation due to 


surface roughness. 


G.. REFLECTION POLARIZATION 


For thermal radiation at a surface 


0 (8) +a(8) +1 (8) =1 ¢ (6)=a(6) (lus) 


As stated previously, the total amount of energy at a 
surface consists completely of the reflected, absorbed, and 
transmitted fractions. For a perfectly opaque body, 


therefore 


ano) (Oo) 1 (Ey 


The emissivity (€) was treated in the previous section; the 


reflectance (p) of the sky, sea surface, and target must now 


be considered. Absorption will have an effect on the 
reflectance for each surface considered and will act to 


degrade the reflectance. 


26 


1. Sea Reflection Polarization 

The application of the Fresnel Equations to a Cox-Munk 
sea surface slope distribution can be shown to predict 
polarization in reflectance as well as emittance with the 


f£OLlOowing melationsie 1]. 


ope emus 120) 


ee A ae) 
GIS). 


P= 1/2 (P+) 


g=1.9 


The sea surface acts as an opaque body. Values of € may be 
obtained as discussed in the previous section. The 
percentage polarization seen due to reflection is defined as 


Till 


= EG 
P= 100-—+—L haaes 


a Sky Reflectance 
As stated in the emission section, the sky acts as a 
poor emitter near the zenith angle due to the low optical 


depth. Because of the high content of atmospheric water 


Zl 


vapor and aerosols, it acts as an excellent absorber. By 
equation (13) it 1s therefore a poor reflector. 

i Target Reflectance 

The characteristics of target reflectance follow the 
same behavior as the wave facet reflectance. The difference 
in the magnitude of the polarized components is a function 
of the index of refraction of the surface material of the 
target. As seen in previous cases, maximum polarization 
effects due to emission, and thus polarization will occur 
nearer to grazing incidence angles. The net polarization 
will be dependent on the balance between thermal emission 
and reflection of incident radiance, related through Eq. 
(15). For the purpose of the experiments performed to 
demonstrate the different polarizing imaging techniques, 
target distance and lack of detailed geometric orientation 
of surfaces precludes exact measurement of target surface 


polarization; however, qualitative effects will be shown. 
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IV. DATA ACQUISITION TECHNIQUES 


To analyze the target radiance as it is affected by 
temperaturer, polarization, and other factors, it is 
necessary to employ thermal imaging equipment for recording. 
The split-field polarimeter and internal filters use the 
AGA-780 thermal imaging system for data acquisition. In 
the case of the internal filters which were designed to fit 
in only the AGA-780 Thermovision, the choice is necessary, 
but the split-field is made to adapt to a number of 
different sensors, particularly the AMBER 4128 focal plane 
array system. All the experimental data taken to be 
analyzed was taken with the AGA-780 Thermovision. 
Accompanying software necessary to digitize and display 
radiance data consisted of the AGA CATS software and 
Interactive Data Language (IDL) programs. The following is 
a brief discussion of the data collection devices. 

A. AGA-780 THERMOVISION 
The AGA-708 Thermovision consists of a two channel 


infrared scanning detector system. The AGA can operate in 
both the 2-5.6um and 8-12um bands using liquid nitrogen 


cooled InSb and HgCdTe detectors, respectively. Scanning 1s 
accomplished by two rotating prisms which scan the image 


across the detector, yielding 1.1 mrad geometric resolution. 
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The input lenses for the Long-wave infrared (LWIR) used for 
the experiments consist of anti-reflection coated Germanium. 
The LWIR band was used for all imaging. The optical system 
is £/1.87, with a 7x7 degree field of view. This f number 
was a key design input for the split system to match the 
aperture of the AGA. The AGA-780 is mounted on a common 
plate with the split-field system with a sliding mount to 
facilitate easy changing between the LWIR and MWIR band 
detectors. Figure 11 in chap. 7 shows the arrangement used 
with the split system. 
Ere AGA CATS SOFTWARE PACKAGE 

The analog data provided by the scanning detector is 
digitized for display and measurement by the AGA CATS 
software. The CATS software uses an 8 bit A/D converter to 
convert detector intensity into a 16 color display for 
real-time viewing on a computer screen. Many different 
functions are available to the user such as storing frozen 
image data, image manipulation, temperature readout and 
emmissivity compensation. When using the CATS program to 
provide temperature readings of a target, it iS necessary to 
calibrate the system to match target radiance with 
temperature. This is accomplished by recording the thermal 
level seen by the A/D conversion and fitting a curve to the 


data through a least-squares method. The curve is 
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Peo, and C are determined by Ehe Gurve £lLtting. = Ose1s. che 
offset correction up or down, and T is the temperature 
(ordinate) of the plot. Thermal Level value range 
determines the abscissa. Plots for system calibration for 
the internal filters and with the split-image system added 
are shown in Appendix D. 

The most important feature of CATS is the ability to 
store frozen image data for later analysis. To make 
detailed measurements and perform data manipulation, further 
software is reguired with greater flexibility. 

c. INTERACTIVE DATA LANGUAGE (IDL) SOFTWARE 

In order to process the image data to a level 
comensurate with the accuracy in measurement desired, a 
software package is required that provides flexibility, 
eegrial resolution, and data manipulation capability. IDL 
allows the user all of these functions and is required to 
give the data processing power needed. It is produced by 
Research Systems Incorporated, Boulder, CO. Data that is 
stored by the CATS program is read into IDL in array format 


where numerical analysis is performed with the data 
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displayed in an image format. Routines are written to 
perform a number of functions and will be discussed in the 


analysis chapter. 
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V. POLARIZATION FILTERING TECHNIQUES 


Use of polarization signature and background features 
for target discrimination implies the employment of 
polarization filters. The technology for large aperture 
polarizers in the IR is relatively new. Most methods 
commonly used in the visible range are not available in the 
DR. 

A. POLARIZING FILTERS 

Polarizing filters have traditionally been constructed 
uSing anisotropic crystals or wire grids. Of these methods, 
the wire grid technique has been applied in recent years 
moing lithographic techniques. THs consTsts™0f a 
conducting grid deposited on a transparent substrate. 

Placed in the optical path of the incoming radiation, the 
filter will pass the component of the radiation with the 
electric field orthogonal to the orientation of the grid. 
The grid behaves like a dielectric sheet in this case. The 
component of the radiation with the electric vector parallel 
to the wire grid induces currents in the wires and is 
reflected. Fig. 9 shows a schematic of this process. A 
restrictions on the grid is that the line spacing must be of 
the order of the wavelength of the incident radiation 


(Ref.6). Such filters for the infrared are available with 
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clear apertures up to 45mm. Typical transparent substrates 


in the infra-red are materials such as Germanium, Silicon, 


and KRS-5. 
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The location of the filter in the optical path has a 
large impact on image attributes. If a flat reflecting 
surface is placed in front of the collimated lens, it is 
possible for the cold detector to 'see' itself in the image. 
This is referred to as the 'Narcissus Effect' (Ref.23),; 
which perturbs the image greatly. Tilting the polarizer may 
alleviate this problem but will reflect radiation into the 
image from outside the field of view. Another drawback to 
using an external filter is the size required to cover the 
aperture defined by the input lenses. It is necessary to 


cover the entire aperture which is a function of the system 
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f number and focal length. Benefits of using an external 
filter are ease of insertion/removal and grid orientation 
changing. 

The large external filter was used in the MAPTIP exper- 
iment in October 1993. The large external polarizer pro- 
vided excellent data for the determination of the degree of 
polarization of the targets and backgrounds and of contrast 
improvement (Ref.1). Conclusions showed percent polarization 
of vertical sea-water emission up to 83% in the long-wave IR 
band. Contrast improvement from emission polarization 
filtering was as high as 36.3%. Large improvements such as 
this justify the utility of the external polarizers and show 
their effectiveness. 

Dr. Internal Filters 

Filters can also be placed inside the AGA-780 imaging 
system between the main input lenses and the detector. This 
allows hands-off operation and smaller filter size. In 
systems which are configured to use filters in the optical 
path for other types of filtering, small polarizing filters 
can be substituted. These filters will typically be de- 
Signed into non-imaging locations. ‘'Narcissus' spots, 
filter mounting and alignment problems, and other complica- 
tions can thus be avoided with this method. Implementation 


can be accomplished with minor system modification and cost. 
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Installation in an existing filter wheel requires a polar- 
izer for each polarizer orientation to be used, adding to 
system expense. Internal filtering is a highly desirable 
option and experimental use will be discussed in more detail 
later. 

An alternative method of internal polarizing applicable 
to focal plane array systems is placing an independent 
filter closely against each element of the focal plane 
array[6]. Each element of the array can then have any 
preselected orientation of linear polarization built into 
the focal plane mask. An added benefit over rotating 
polarizers is the simultaneous recording of all the selected 
polarizations without manipulating a filter selector. 
Disadvantages are the pixel offset of the image for each 


polarization, and the inherent loss of resolution. The 
practicality of this technique for the 3-5um band has been 


demonstrated in a NASA-Supported program (Ref. 6). 
B. SPLIT FIELD OPTICS 

A method of polarization in which the polarizer is 
internal to an optical system but external to the imaging 
system itself is the technique of Split-Field polarization 
filtering. The concept of the split field system is that 
orthogonally polarized images can be extracted from the 


field-of-view simultaneously with one polarizing filter and 
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presented in adjacent sections of the final image field of 
view. This is accomplished by splitting the beam with a 
polarizer. The Environmental Research Institute of Michigan 
(ERIM) has developed such a system which uses Brewster angle 
reflection to separate one plane polarized component from 
the beam, extracting the orthogonal component from the 
transmitted component through the Brewster plate (Ref. 5). 
Since the Brewster plate is tilted, a reference is needed to 
keep uncalibrated background from contaminating the image, 
as in the NPS system. Commonality exists between the ERIM 
system and the NPS system. However, the NPS system defines 
the apertures of the two channels and brings them to 
adjacent halves of the final image plane. The ERIM system 
focusses the two channels onto separated areas on a single 
integrated focal plane array. 

The Naval Postgraduate School, Naval Academic Center 
for Infrared Technology (NACIT) has constructed and begun 
testing a split-field system. Polarization is attained 
through a wire grid polarizer tilted to the incoming 
radiation with polarization orientation obstained from grid 
transmission and grid reflection. When considering options 
for selecting an imaging system for the split-field optics, 


the scanning AGA was chosen. The only focal plane array 


available for use at this time operates in the 3-5 um band. 
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Because of the nature of the targets and environments where 


the split system will be used most, it was decided to first 
build the system to operate in the 8-12 um region. The 


existing hardware and familiarity using the AGA coupled with 
IDL software and established data analysis techniques made 
it a good choice for examining the new system. The trade to 
go to the AGA is in reduced frame sampling rate and 
resolution. A detailed explanation of the split-field 
technique employed, a system description, and polarizing 


characteristics of the system will follow later. 
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VI. INTERNAL POLARIZING FILTERING WITH THE AGA-780 


Internal filtering of the radiation after it passes 
through the imaging system input lenses can be accomplished 
with the AGA-780 Thermovision Scanning Imaging system by 
using an existing filter wheel. The original purpose of the 
filter wheel was to hold various accessory spectral filters 
in the optical path before the radiation reached the 
detector. Fig. 10 shows the internal layout of the AGA-780. 
One of seven positions of the wheel is selected by rotating 


the knob on the front face of the unit. Care must be taken 
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Figure 10. Internal Layout of the AGA 
Thermovision (model with only one IR Band). 
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to ensure that no internal circuitry compensates for an 
assumed filter present at the position of the polarizer. 
Polarization measurements with this system require 
sequential images with alternate filter wheel positions 
inserting vertically and horizontally polarized filters. 
This necessarily involves time delays. 
A. ALUMINUM/KRS-5 FILTERS 

The filters used consisted of an aluminum grid 
Superimposed on a 9.5mm diameter x 3mm thick round KRS-5 
substrate manufactured by Graseby-Specac; Suffolk, England. 
The performance specifications provided with the filters are 
nearly identical for both and will be treated as such for 
explanation. Appendix A shows the performance curves for 
the two filters. The various traces numbered 1-3 represent 
power transmitted with the grid vertical, horizontal, and 
with the grids from two filters crossed, relative to the 
electric vector of the incoming radiation. The determina- 
tion of the perpendicular and parallel transmittances for 
the filter involves use of a spectrometer which itself 
includes internal polarizing reflections. Thus it must be 
accomplished by measuring the transmittance twice, with the 
grid in two orthogonal orientations. These two measured 
transmittances are designated E, and E,. The corresponding 


instrument transmittance is measured without the grid in 
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place as E,, and the transmittance for two crossed grids as 
E,. From these quantities the grid transmittances are 


derived as 








These quantities are specified by the manufacturer for each 
filter. The expressions are approximate values derived by 
Mmewngd,. CT dla (Ret . aS). 
By FILTER INSTALLATION IN THE AGA-780 

Figure 10 shows the internal layout of the AGA-780. To 
install the filters, the appropriate holes in the wheel were 
enlarged approximately 0.5mm in diameter. Two filters were 
then installed in filter positions 4 and 5, carefully 
oriented with passage axis vertical and horizontal, 
respectively. The filter wheel was then installed in 
accordance with the manual. Extreme care was taken to 
select positions in the filter wheel which have no 
electronic filter compensation in the system. The manual 
denotes which positions these are. Once the system 
modification was complete, new calibration constants were 
computed and installed in the CATS system. Calibration 


constants are listed in Appendix D. 
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VII. THE NPS-NACIT SPLIT-FIELD POLARIMETER 


The split-field system approach possesses many 
advantages over external and internal filter placement 
mormal to the path of bhe zadiation, particularly im 
elimination of time lags and registration problems between 
orthogonal polarizations. The Naval Academic Center for 
Infrared Technology (NACIT) at the Naval Postgraduate School 
(NPS) has constructed and tested a split-field front end 
optical system for input to thermal imaging systems, 
specifically matched to the AGA-780 and the AMBER 
Engineering-4128. The Split-Field Polarimeter is the result 
of a conceptual design by Professor Alfred W. Cooper. The 
actual construction and design optimization was performed by 
Professor Eugene C. Crittenden and LT Peter H. Heisey, USN, 
as assistant. The general features of the NPS-NACIT Split- 
Field Polarimeter are: 

1. Simultaneous display and recording of orthogonally 

polarized images. 

2. Utilization of entire field-of-view for display. 


3. Adaptability to different imaging systems. 


4. Dual-band (8-12um and 3-5um) operation using 
interchangeable filters. 


S.  f/71287, 99nm f¢6¢a | SlenGeheorEtcs . 
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As mentioned Wn Chapt. iV for other systems with split-field 
optics, the 'Narcissus spot' is avoided by having the filter 
tilted relative to the path of the radiation. A background 
plate is also needed to give a uniform reflected radiation 
source into the images. The specifics of the system and its 
components and design will be discussed. 
A. SYSTEM DESIGN 

The split-field optical system consists of input and 
output compound lenses, polarizing splitter, and front- 
surface aperture limiting mirrors. The components are 
arranged in a vertical plane with respect to the ground. 
The system takes the collimated input radiation, focuses it 
down as it is split into two orthogonally polarized 
components, and folds and recombines the two split halves of 
the image while recollimating the light for processing by 
the thermal imaging system. In this process, the optical 
aperture is redefined so that only the scan lines of the 
central half of of each image are retained. The remaining 
3.5 x 7 degree halves are combined in adjacent halves of an 
intermediate 7 x 7 degree image plane. The 7 x 7 degree 
field of view of the optical system is thus split into two 
channels with orthogonal polarizations. Each channel is 
aperture limited into a 7x3.5 degree field to compose half 


of the final image. The original radiant flux in the image 
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1s preserved in this way. An overall picture of the system 
in use with the AGA-780 Thermal Imaging System is shown in 
Fig. 11. In this setup, the entire imaging system is 
mounted on one lightweight mounting plate. Fig. 12 shows 
the inside of the box with all components exposed. Aluminum 
was used for all mountings and plates because of low weight, 


corrosion resistance, and easy machining. 
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ProGaror aaa, Top View of AGA with Split-System. 


45 


wae was gy 


at 








Figure 12. Side View of AGA with Set eee les 


a 1 os tow 
RREKOKKAnMRe BAKK 
ME RAENS ER CATER ROG 
ARMM MAKR ORE AY 
rem es Ones ere Mary 
Ae AE RAS OR ERERE Be KY 6 BAL 
AKT AAS AY CU REKAS > 
RARKUAAYXEKAMMKK DX 
Pa NEY KNROR GON OKAM Ge 
WH Coe  KRUKHHKKOKRAM VO KE 
£RY Koo RV URRR OMAK 
4, Ween cK ZNUES 
KARR ex 
MOM KAMAKEKK 


otk wr) 


xa 
ox 





Figure 13. Front View of Split-Field input 
lenses with AGA in background. 
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availability, corrosion resistance, and easy machining. The 
front view in Fig. 13 shows the exposed input lens. The 
long wave input lens for the AGA system is hidden by the 
backside of the system. The short wave lens not being used 
by the AGA is to the side and may be slid into place behind 
the optics when a short-wave polarizer is installed. Figure 
14 shows the image orientation and focal plane within the 
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Figure 14. Schematic of Ray Paths in Split Optical System 


The incoming rays encounter two Germanium lenses (labeled 


lens system) and are refracted to a focal plane. In the 
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bi 0 CeSiS msec wore 7ainiGmic Diem SO kits  thesbeameinto two 
paths with one beam vertically polarized with respect to the 
ground (transmitted beam) and the other reflected with 
horizontal polarization. The transmitted beam is reflected 
to a focal plane at the edge of a small mirror used by the 
reflected beam to be turned in the same direction (Fig. 14). 
At this point, both beams have been recombined in a single 
focal plane. The beams then diverge and are reflected by the 
last mirror into two lenses identical to the first lens 
system (labeled lens system). The radiation emerges 
afocally so that the AGA imager, positioned immediately 
behind the system can process the image for display. 
Bi. SPLIT-FIELD OPTICAL COMPONENTS 

The split-field system consists of 3 major components 
as stated in the previous chapter. The specific 
characteristics and operation of these components will now 
be addressed. 

- Lenses 

The lenses used by the system to generate the split 
image for input to the thermal imaging system consist of two 
identical assemblies each containing two germanium anti- 
reflection coated lenses. One assembly is used for the input 
and the other for the output side. The output assembly is 


reversed, end for end, relative to the input assembly. 
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The inner surface of the first lens of the input assembly 
and the inner surface of the last lens of the output 
assembly are aspheric. This reduces aberrations and provides 
a minimum spot size of 1.1 mrad, within the resolution of 
the AGA imager. The lens surfaces are coated for 
antireflection in both the 4um and 12um wave bands. This is 
necessary because the high index of refraction of germanium 
(4.0) would produce large reflection losses and scattered 
radiation if the surfaces were left uncoated. The overall 
focal length of each of the two lens assemblies is 99.2mm 
with an f£ number of 1.87 

an Polarizing Splitter 

The most important component of the split system from 
the polarizing standpoint is the Al/Ge polarizing splitter. 
The custom designed splitter consists of a germanium 


substrate with an aluminum grid deposited on one side. The 


grid consists of .25um wide aluminum strips separated by 


~-40um for an overall period of .65um. The substrate is a 
round disc, 50mm in diameter x 2mm thick. Both sides of the 
splitter areanti-reflection coated for the 8-l2um band. It 


1s necessary to replace the splitter for the short-wave IR 
band, since the vendor found it impossible to coat the 
element for both bands. The same analysis done for the KRS- 


5 polarizers in Chapt. V was performed for the Ge splitter. 
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Data provided by the manufacturer (Graseby-Specac) is valid 
for normal incidence of the optical beam on the splitter. 
Values for K,, K,, and P are .81, 0.0, and 1.0, respectively 
at 10pm. This means the splitter passes 81% of the linearly 
polarized light perpendicular to the grid orientation and is 
~100% polarized. The manufacturers performance curves are 
provided in Appendix A. Extensive analysis was done to 
verify operation of the splitter; this will be addressed in 
section C. Results showed that there is a wavelength 


dependency of the radiation fraction that is reflected and 


transmitted. It is found that wavelengths between 7um and 


9um favor reflection and the opposite is true above 9um. 
This is due to inefficiencies in the polarizing grid and 
anti-reflection coating. Results obtained by system 
calibration show this and will be discussed later. 

3. Gold Mirrors 

The folding operation of the optical system is 
performed by four gold first-surface 'protected' mirrors. 
The mirrors were purchased from Edmund Scientific Co. and 
are advertised to be 1/10 A flat at .5yum. Gold was chosen 
because of its high reflectivity, and thus, low emissivity. 


Figure 15 shows a plot of the emissivities at 9.92 um for 


gold at a range of incidence angles. 
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As expected, significant polarization effects occur only at 
angles close to grazing incidence. Most reflecting angles 
of the mirrors are near 45 degrees, so these effects are 
small (<1%), but must be considered. 
c, SYSTEM CHARACTERISTICS 

The combined effects of the presence of the components 
in the optical path and mechanisms by which they perform 
their functions must be quantified to predict overall system 


performance. Irregularities in refracting and reflecting 
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surfaces, polarized reflection and emission, and system 
background radiance are effects which influence the measured 
radiance and resolution of the target. Experiments have 
been performed to characterize these effects to assist in 
predicting final image quality. 

Ly: Resolution 

To obtain the full resolution possible by the AGA-780 
when used with the split-field system, effects induced by 
imperfections in the optics of the split-system must be 
minimized. Upon delivery of the splitter, visual resolution 
was tested using a standard bar chart. Reflection of the 
bar chart off the surfaces of the splitter showed no loss in 
resolution from direct viewing of the chart. It follows 
that imperfections in the surfaces would induce less effect 
in the longer wavelengths and would not be considered. 

A test was also done on the combined AGA-split system 
to discern if resolution was degraded by the presence of the 
split-field system. A sharp edged (cold) plate was placed 
in front of a black-body thermal source. Using the AGA 
optics alone, measurements were taken of the number of 
pixels in the image that were required to resolve the sharp 
edge between hot and cold surfaces. Along the direction of 
image raster scan, the resolution is about 2-3 pixels. 


Perpemdicular to the raster scan the resolution is slightly 
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better. Parallel to the raster scan the resolution is 
apparently limited by the frequency band pass of the 
amplifier electronics and the response time of the detector. 
In the direction perpendicular to the raster scan, the 
resolution is slightly higher because the response times of 
the detector and amplifier are not involved. In both 
parallel and perpendicular scan directions, no loss of 
resolution occurred in the image when the split-field optics 
were added. 

2. MTF 

Modulation transfer function tests were performed for 
the AGA system with and without the split-field. The cutoff 
frequency for the AGA in the direction of raster scan is .3- 
.5 cycles/mrad which gives a resolution of 2-3mrad. The AGA 
system manual gives a system resolution of 1.1lmrad. This is 
apparently the geometrical resolution of the AGA optical 
system. The MTF data show that the overall system has a 2-3 
mrad resolution. There was no loss in spatial resolution 
with the split-field system; 

aS System Polarization Performance 

The main driver in the performance of the system is the 
polarizing splitter. Establishing numerical values for 
transmissivity, reflectivity, and absorption of the splitter 


for polarized and unpolarized target radiation requires 
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measurement of all radiance sources throughout the image and 
split-field system. 

a. Spectral Polarizer Performance 

The manufacturer of the polarizer provided 
spectral transmittance data for the polarizer at normal 
incidence to a reference beam. For use in the split system, 
it 1s necessary to know the spectral transmittance and 
reflection for 45 degree incidence. 

To measure the transmittance, the polarizer was placed 
at 45 degrees to the incoming beam in a Perkin-Elmer Model 
337 Infrared Grating Spectrophotometer. The photometer 
measures the characteristics in the same manner as the 
spectrophotometer used by the manufacturer. Because of 
polarization induced in the beam by the grating, the method 
described by Young et al. (Ref 13), is used to compute the 
transmittance. This is the approach outlined in Chapt. VI. 
Appendix B shows the plots generated by the polarizer at 
normal and 45 degree incidence. Agreement with the 
manufacturers values are within 3% for normal incidence. 
For wavelengths less than 10.5um, the transmissivity is 
higher than at normal incidence. This is probably due to 
the fact that the effective grid spacing is reduced by 
cos(45), when the grid is perpendicular to the plane of 


incidence of the beam, thus moving the curve left to shorter 
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wavelengths. It is necessary to mention that placing the 
filter in such a manner as to reflect ambient radiation into 
the detector does not contaminate the readings. The 
intensity radiated by a typical room temperature object 
compared to that of the IR source in the device is on the 
order of 0.001. 

Obtaining performance for reflection polarization 
required altering the conventional configuration of the 
polarizer to the beam to utilize the front surface as the 
effective transmitter. Figure 16 shows the arrangement used 


to obtain this. 
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Data for the linear polarization reflection efficiency is 
provided in Appendix B. Reference 13 provides a derivation 
of the measurement technique to obtain the transmittance of 
a polarizing grid. The difference in finding the efficiency 
here is that reflection is substituted for transmission in 
the equations. The radiant intensity with the polarizer 
removed was obtained by placing an aluminum mirror in place 
of the polarizer and resetting the 100% transmission 
indicator. The performance plot in appendix B shows this 
trace across the top of the graph. The designation for the 
reflected efficiency is K,,. The equations are analogous to 
those described in chapter VI. It was not possible to 
measure K,, (orthogonal polarization reflection) because of 
the unique arrangement and nonavailability of another 
Similar polarizer. 

The degree of reflected polarization perpendicular 
to the grid, can be described as follows. Any radiation 
impinging on the grid with E vector perpendicular that is 
not transmitted must be reflected due to fact that the grid 
width is small compared to the wavelength, precluding 
current generation. This reflected component is the 
complement to the transmitted fraction of the radiation. 


This theory will be demonstrated in the following section. 
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De Polarizer Performance Utilizing Radiance 
Analysis 


(1) Theoretical Representation 

An alternative method for meaSuring the fraction 
of polarized radiation transmitted and reflected by the 
polarizing splitter is to take radiance readings of the 
radiation passed with a known source polarization 
orientation. However, complications arise with this method 
when target temperature is near the temperature of the 
optical system, because the background irradiance that 
enters the detector is of the order of the target 
irradiance. To discriminate the target irradiance from the 
background, a mathematical description of the background 
effect is required. 

For the case of unpolarized radiation incident on 
the optical system, the following equations apply: 
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I and I, represent the components of target irradiance 


perpendicular and parallel to plane of incidence of the 
polarizing splitter. These are what is actually measured by 
the thermal imaging system. For the entire analysis the 


convention used is that the polarizing splitter grid lines 
are perpendicular to the plane of incidence. The €,,, terms 


represent the radiant exitance of all the components of the 
system when the incident intensity into the system ~0. This 
is simulated by directing the optics at a pool of liquid 
nitrogen. The subscripts T and R designate the transmitted 
and reflected channel through the split-field optical 
system. T denotes the radiance fraction which passes 
through the polarizing splitter and is displayed as the top 
half of the split image; R denotes the other channel 
reflected by the splitter. The perpendicular and parallel 


subscripts describe the orientation of the E vector for the 
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irradiance relative to the plane of incidence of the 
splitter. For sexamplepeunpolari zed, radiation sseeommenly 


represented by 


meme 9/62) bate (lin 5/ 2) (20) 


total ( total ats total 


To eliminate unknown terms in Eq. (19), measuring 
a purely linearly polarized source in two orthogonal 
orientations yields a simpler set of equations. For the 
first case of purely polarized target irradiance 
perpendicular to the plane of incidence (perpendicular to 


the grid), 


I,~f, ay ane (21) 
r,t, a re, te, 


Ry 


Mime terms denote Ehat fraction which gocs intogene 
respective channels from splitter reflection and 
transmission performance. The target consisted of a 
memarizing fidter in Eront of a black beady source; 
therefore, there is some filter emittance that must be 


considered. The added €, terms represent the emittance of 


the polarizing source for each channel, for each orientation 


oo 


of the target producing filter. It can be seen that there is 
a fractional contribution to each channel to be measured. 7T 
and R. This is because it is known that the transmittance 
of the grid is not 1.0. For the case where the target 
radiance is plane polarized perpendicular to the plane of 


incidence (parallel to the grid), 


I_=e_te 
Te i E 
B Ds (22) 
fae Ire, te r, 


a 1 ue 


There are few terms seen in the T channel because it is 
known from the spectral performance curves that the 
transmitted E components are ~100% polarized. Equation 21 
shows the intensities seen in the two channels with a 
vertically, parallel (p) linearly polarized source and eq. 
22 with 90 degree rotation of the target E vector, s 
polarized. 

(2) Experimental Radiance Measurement 

To compute the performance of the system using the 
equations in the previous section, it is necessary to take 
careful radiance measurements of each term in the equations 
with a linearly polarized target. Appendix C contains the 


values obtained in the laboratory for these measurements. 
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To compare the performance of the polarizing splitter 
channels, the values in Appendix C are plugged into egns. 
(21) and (22) and then subtract the equations from each 


other. The relationships obtained are 
(f,-£,)1,—21.9 


a eae 


Measurements showed that there is no dependency of target 
intensity on polarization orientation, so IL-l, . 


Therefore, 


are) 


—_+—__1-=0.98 
F (23) 


R, 


This equation allows calculation of the reflected fraction 
of the incident radiation which is normally passed through 
meeche grid. Whe implications Of Eq. (23) are that the 
transmittance of linearly polarized irradiance perpendicular 
to the grid is always greater than the reflection efficiency 
of the orthogonal component from the grid front. This is 
seen with the performance curves in appendix B. It also 
allows correlation of splitter performance using the 
radiance method with that seen in the spectral performance 


curves. To compare the radiance method, it will be 


ol 


assumed that the non-transmitted E vector perpendicular to 
the grid undergoes Snell reflection and f,+f,=1. The 
spectrophotometer data measured in the previous section 
shows that at 9.52um, the above expression = 0.97. 

Agreement is seen between the two methods at that 
wavelength, but the ratio varies strongly at other 
wavelengths. Schematic representation of what happens at 
the splitter within the optical system is shown in Figs. 1/7- 
19. Figure 17 indicates the fractional performance in the 
LWIR for unpolarized input. For this case, 58% of the 


background is reflected into the transmitted (T) channel and 


Reflected component 
ONPOLARTZED INCIDENT RADIATION 0.42 Ib 


Itgt  ws----+- an Polarizer 
0.58 Th 0.41 Itgt 
pe 59 Itgt l 
see Cae Transmitted Component 
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Figure 17. Schematic of polarizing splitter 
for unpolarized incident radiation. Grid 
perpendicular to plane of incidence. 
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combined with 41% of the target irradiance. The reflected 
(R) channel combines 42% of the background with 59% of the 


target irradiance. Figure 18 indicates the fractional 


Reflected component 


HORIZONTALLY POLARIZED INCIDENT RADIATICR 


Itgt 


@ 





Transmitted Component 


Background Plate 


Figure 18. Schematic of polarizing splitter 
tom horazontally polarized mmerdent 
radiation. Grid perpendicular to plane of 
incidence. 


performance in the LWIR for horizontally polarized input. 
For this case the transmitted beam contains 58% of the 
lrradiance from the background plate and essentially no 
target radiation. Figure 19 shows the equivalent situation 
for vertically polarized input, which shows 83% of the 
target irradiance with 58% of the background irradiance in 


the transmitted (T) beam, and 17% of the target irradiance 
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Reflected component 
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Figure 19... Scnememitewe enpelarizing splierer 
for vertically polarized incident radiation. 
Grid perpendicular to plane of incidence. 


with 42% of the background plate irradiance in the reflected 
(R) beam. 

The radiance analysis of the polarizing splitter 
performance provides a gross analysis and verification of 
the data provided by the manufacturer and by the grating 
spectrophotometer. This representation of transmittted and 
reflected target irradiance must be considered when 


processing split-field images. 
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VIII. IMAGE ANALYSIS 


To prove the concepts presented and to demonstrate the 
etfects of polarized fisjtering on imagese@ontrast, 22 1c 
necessary to perform a preliminary image analysis on images 
recorded with the thermal imaging system in two 
configurations; internal polarizing and split-image optical 
system polarizing. Two separate experiments were conducted 
to demonstrate system performance. 

In March 1995, 1338 images were taken during a field 
experiment with the R/V POINT SUR from the Monterey Bay 
Aquarium Research Institute marine operations building at 
Moss Landing, CA. All of the polarization filtering was 
done internally. In August of the same year, images were 
recorded using the NACIT Split-field Polarimeter from the 
Naval Postgraduate School beach property of the USS John 
McCain, an Arleigh Burke Class DDG. The images recorded 
show contrast improvement effects from polarization 
imo tering. 


A. DEFINITION OF CONTRAST AND TARGET-BACKGROUND RADIANCE 
DIFFERENCE 


When an observer views an image, the critical parameter 
for detecting and recognizing the target amidst the 


background is target contrast. Numerous definitions are 
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currently used to definewcontrast. Common examples of these 


are [23] 
Contrast= IV argeee pact erenne 2 Contrast= EY aera ee ore and 
a ni, <a ae 
(24) 
background -. N 
target background 


Here N represents radiance. It is desirable to compare the 
contrast of images which may be separated in time. In this 
case, changes in system gain cause the apparent digitized 
image values to change for similar viewing conditions. 

Also, changing instrument settings changes N for every image 
when the controls are altered. Therefore, when target and 
background conditions vary little between images, the 
difference in N may be a more accurate indicator of 
contrast. Differences in filter transmittance cause the 
conventional contrast equations to be in error so that any 
artificial thermal level changes induced by this are ignored 
for filters with the same calibration curve slope. This is 
the case for the internal filters (Appendix E). To 
accommodate this, the analysis will only consider 
differences in the apparent radiances of the target and 
background. The contrast parameter to be measured is termed 
Target-Background Radiance Difference (TBRD). This is the 


numerator of both contrast equations in Eq. (24). 
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fBRD,— AiG Swe > ~~ <N > (25) 


target background 


The average values of N are used because of the varying 
composition of the image with target distance, system noise, 
and measurement points in the image. If the image radiation 
1s polarization filtered, the resultant TBRD is a function 
of the amount of filtering only if all optical components 
and scene characteristics remain constant. It has been 
shown that many different components of an image can display 
different degrees of polarization. If background irradiance 
1s suitably filtered more than target irradiance, TBRD will 
increase. 

When comparing two images, it is useful to quantify the 
increase in image contrast improvement. The TBRD 


improvement factor is defined by 


AN 


TBRD. =AN_ =—— 
improvement imp A N, 





B. FIELD MEASUREMENTS WITH INTERNAL FILTERS 
From May 16-23, 1995, polarrzea@sthenmal imaging 
experiments were conducted from the Monterey Bay Aquarium 


Research Institute using the R/V POINT SUR as target. 
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1338 images were recorded using internally filtered and 
unfiltered images in the 8-12um wave band. The ship 


Preseobee Salish ecetseWaetuticcges Varyingetronm.25enmi to 
the range limit of the day. Appendix E carries a complete 
listing of the images and a small picture of each with the 
filename shown on the bottom of each image. Vertically and 
horizontally polarized images will be compared for target- 
background radiance difference. 

cs Visual Image Comparison 

Selected pairs of images demonstrate the visual TBRD 
and contrast improvement between vertically and horizontally 
polarized images. Unpolarized images were not considered 
for this comparison because of the nature of the data 
analysis process, and the change in system settings 
necessary to obtain useable images due to the changing of 
system gain and altering the byte level of the images 
Significantly. The definition of TBRD would not be valid 
for comparison to the polarized images and further 
compensation to TBRD would be necessary. The images are 
selected from the compilation listing in Appendix E. The 
contrast between the two orthogonally polarized members of 
each pair of images in Fig. 20a and 20b shows the lower 
apparent radiance of the sea water background in the 


horizontal polarization. It is necessary to mention that 
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the images were printed on a Hewlett-Packard LaserJet 4M 
600dpi printer. Differences in the images are more easily 
seen in a grayscale color table on an SVGA monitor. The 
bright area seen in the top right corner of each image is a 
constant temperature source created by inserting a mirror in 
the field of view, reflecting a calibrated black body 
source. This provides an unpolarized reference in each 
image to compare filter transmittance and provide 
temperature calibration confirmation. The source is not 


discussed in this analysis. 


ow 
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Figqune. 7 UaweG@empaGisicn Of Venmealiyeand horizontally 
polarized LWIR Band (8-12ym) images of RV POINT SUR. 
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Paoure 20s Comparison OL verumealy ane morlzonvaldy 
polarized LWIR Band (8-12pm) images of R/V POINT SUR 


ae Numerical TBRD Measurement. 

The images presented in Fig. 20 were extracted groups 
of images with the same aspects, ranges, and path 
conditions. The average radiances are calculated by 
averaging over the target and background areas within an 
area around the ship. Appendix F contains the IDL code to 
perform this operation. To compare the TBRD for a number of 
images, the images must be taken with target and background 
conditions as constant as possible for trends in the data to 
appear. Figures 21-23 show TBRD measurements made on 
different sets of images grouped by ship range and aspect. 
The vertical axis gives the TBRD in dimensionless units. 

The 8 bit A/D converter sets the maximum thermal level 
reading based on system gain to 255. This irradiance 
reading is independent of system calibration settings and 
measures the raw data from the detector with no compensation 
for temperature accuracy. The horizontal axis gives range 
computed from ship image height measured in pixels. Smaller 
targets appear to the right side of the plot, which shows 
TBRD as a function of range. The general trend of TBRD 
reduction with target size is due to increased path 
attenuation and path radiance with range. The legend on the 
plots denotes the polarization and whether TBRD is 


calculated above or below the horizon. The average TBRD 
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Figure 21. TBRD Improvement-Starboard Aspect 
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Figure 22. TBRD Improvement-Port and Stb Aspect 


74 


dZI§ diyg Wor ‘SOTA, [BOWNBNY ‘asuvy diys 


(9° SS" QS" SP" OV So" 
cc © © + f 
o ? + 
2 a zi 
? 
0 O O CJ 
O) O ‘ 
* x > : 


CF T=<Mowesoidun QQ L> 
A=90 H= +e 


% + OSLIDY “0[g 


Os woxoq eacgy 


ISURY ‘SA CUAL 0L09T-0S09TS TN 


Od 


OF 


09 


08 


001 


0c 


JDUaIIJJIq BJOUPIPeY PUNOASHOegG—JoAIeL 


Figure 23. TBRD Improvement-Bow Aspect 
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improvemement on each plot is calculated only for the below 
horizon case against the sea water baGKGmeounceur lt. ais 
expected that TBRD against the sky should be nearly the same 
in both polarizations. This data does not contain any 
meteorological, target temperature, or detailed Shay 
position data which is available for analysis improvement. 
Very general trends are seen for TBRD improvement, but 
are not as apparent as the average improvement calculation 
meenGChS In thewplot. “Ship size uvaries greatly in the plome 
due to variation in the numbers of pixels recorded in the 
vertical direction and is most pronounced in the images 
taken in the near ranges. It is seen that there is some 
polarization effect in the above horizon TBRD. This is 
probably due to horizon selecton in the image and possible 
polarization of atmospheric radiance close to the horizon. 


Analysis was not performed on the polarization of the sky. 
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Ce FIELD MEASUREMENTS WITH SPLIT-FIELD POLARIMETER 
Preliminary test and evaluation imaging was done on 3 
August, 1995 with the split-field polarimeter optical 
system. The target is the USS John McCain, an Arleigh Burke 
class destroyer present in Monterey Bay for a port visit at 
the Naval Postgraduate School. Thirteen split-field and 10 
test images were taken of the ship and black-body calibra- 
tion sources. The target was at a stationary range of 
1.85km, determined by target geometrical size. All images 
are port aspect. All of the images taken are listed by 
filename in Appendix G. Because of the similarity in image 
composition, only one image will be analyzed in depth. 
Figure 24 shows the image chosen which is optimized for 
focus, level gain setting, and system adjustment. Images 
are recorded inverted by the AGA and this image has been 
artificially flipped upright for viewing. The top image is 
the reflected component of the polarizing splitter, 
horizontally (s) polarized, and the bottom is the Exrans- 
mitted vertically (p) polarized image. Consistent with 
results seen for the internally polarized images is the 
increase in contrast of the top image. Discussion of 


specific composition of each channel is in Chapter VII. 


Eee) 


Figure 25 shows a vertical profile taken of the image 
through the hot stack. From the known geometry of the ship, 
no significant polarization is expected for the stack which 
is just aft of the large notch in the top of the 
Superstructure. Polarization of emission could exist in the 
ship panels on the bow as they slant outward toward the 
horizontal with height. The bottom of the profile 


corresponds to the bottom of Fig. 24. The x-axis denotes 
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NACIT Split-field polarized image, 
USS John McCain, 3Aug95 
Top-Vertically Polarized TBRD=168 
Bottom-Horizontally Polarized TBRD=2 00 


Figure 24. 
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the 8 bit digitized value of radiant intensities tee etor 
the top image is 200 and 168 for the bottom. This gives a 
TBRD improvement factor = 1.19. This is with reference to 
the average value of the sea water background taken in the 
foreground along the azimuth of the stack in the middle. 

The fog above the image 1s assumed unpolarized. There is an 
overall increase in image intensity in the vertically 
polarized image due to the behavior of the splitter and the 
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Figure 25. Vertical Profile of Split Image, M08030i1 
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The added reflectance from the back of the splitter raises 
the overall level of the vertically polarized image. This 
is another justification for using TBRD to compare image 
contrast differences. 

abe Digital Image Subtraction 

The goal of the split system is to provide simultaneous 
images of the two polarizations which allows digital 
subtraction of the images to visually represent TBRD 
improvement. An IDL program was written to accomplish a 
pixel by pixel overlay of the images of the ships (Appendix 
Ee The program finds the 'hottest' pixel in the image and 
aligns the pixel addresses in the top and bottom while 
Subtracting the values. The resulting image displays a 
value of O (black) where the images had the same intensity. 
The areas in the image where polarization effects occur are 
the brightest. Figure 26 is an image of the top channel 
minus the bottom channel compensated for background 
contribution to the bottom image. This was done by 
adjusting the threshold of intensity for the bottom image. 
This essentially removes any polarization effects seen on 
the ship because of their small effects. 

It is clearly seen that the water shows the largest 


amount of difsference in the image. Irregularities in 
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Figure 26. ‘Digitally Subtracted top and bottom 
split image channels. M0803011, 3Aug95. 


in the image plane show on both sides of the subtracted 


image. Technigues are under development to balance the 


intensity of the two images, remove background 


irregularites, and improve alignment of the two images, 


Od 


but these techniques were not available at the time 
these images were processed. The images are shown here to 


illustrate the general utility and feasibility. 
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IX. DISCUSSION AND CONCLUSIONS 

The split-field polarizing system designed and con- 
structed in this wosk, demenstrates the Utaikity Cliesi museca — 
neous image storing for processing, while preserving image 
resolution. 

Subjective and numerical analysis of the data from 
images taken with the split-field system demonstrated excel- 
lent performance in image quality, contrast improvement, and 
utility. Polarization filtering performance for a sea water 
meaekgqround, for both the split-fielld system and Lor direc 
viewing with interchangeable filters in the AGEMA imager, 
gave results comparable to those obtained in the MAPTIP 
experiments where an external BaFl substrate polarizer was 
used. This was determined by comparing contrast calculated 
in terms of TBRD from previous measurements with values from 
the split-field and internal filtering systems. Both the 
split-field and internally filtered direct viewing systems 
eliminated the narcissus effects that degrade the use of 
external filters. 

Observations of the target scenes show quantitatively 
that: 

1. Digital subtraction of simultaneously recorded split- 


field images is a viable technique for contrast improvement. 


oS 


Leelee soa Peer eC roOn, =verercal pOlariZation emission of 


the sea surface is the predominant component for filtering. 


SemeopliG- field polarized erimkbeerIng Gives preliminary TERD 
improvement of contrast of 1.19. This is driven by the 


splitter performance and background plate effects. 


4. Direct viewing with internal polarization filtering is 
capable of 1.2-1.5 TBRD improvement of contrast for sea- 
water emission, but digital subtraction is not possible 
because of time delay between acquisition of the 


orthogonally polarized images. 


5. Both the split-field and®direct viewing internal f1l]ltegS 
ing techniques provide contrast improvement on the order of 


that seen in the MAPTIP experiments. 


6. The combined system designed and constructed for this 
project provides dual use of the AGA-780 with the split- 
field system and as a direct-view internally-filtered sys- 
tem. This provides two highly effective polarizing systems 
for use in one Heocation with flexibility in system selec 


Cisne 
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X. RECOMMENDATIONS FOR FUTURE WORK 


Further work needs to be performed in all aspects of 
polarized imaging using the split system and internal 
polarizing systems. ~ineMarech, tee an Blectro-eereircal 
Propagation Assessment in the Coastal Environment (EOPACE) 
experiment phase will take place near San Diego CA. In- 
depth demonstration and field implementation of both 
techniques will take place at that time. Specific 
concentration areas for both systems in the intermediate 
time frame include: 

Internal Polarizers: 

-Quantify the spectral transmissivity experimentally 
for each polarizing filter to verify percent polarization 
performance. 

-Re-calibrate the AGA system with the filters installed 


prior to use to give level equalization for both filters. 


Split-image Polarimeter: 

-Conduct laboratory experiments to explore and remedy 
system induced gradient across the field-of-view. 

-Realign Split-field optics to eliminate any 


vignetting, focussing differences between the two channels, 
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and edge effects due to optical surface misalignment. 
-~Install a temperature control system on the background 
plate to maintain constant radiance contribution to the 
image for varying environmental conditions. 
-Insert a small temperature source in the intermediate 
focal plane of the split-field optics to give a constant 


calibration source for real time split channel comparison. 


Overall AGA-780 Thermal Imaging System: 

-Install the new NEC VERSA P/75 75MHz Pentium Laptop 
computer as the display and file storage server to replace 
EhemecunrentaGAhSss SiG sCcomeuber . 

-Install the new 12 bit A/D converter for increasing 
system dynamic range and allowing higher data resolution for 


analysis. 
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APPENDIX A. 


PERFORMANCE CURVES FOR THE POLARIZING FILTERS 


This appendix contains the performance curves for the 
polarizing filters used in the internal and split image 
configuration. The interpretation of the data is included in the 
text of the appropriate chapters. 

The data were provided by the manufacturer of the 
polarizers, Graseby-Specac. The measurements were taken with a 
spectrometer using the actual filter shipped. Where two filters 
were required for the crossed-grid measurement (trace 3), another 


polarizer nearly identical to the one shipped is used. 


Trace 1- transmittance with grid perpendicular to the E vector of 
itj@erdaent radiation. 

feoaee 2— Cransmittance with grid parallel to the E vector of 
incident radiation. 


Trace 3- transmittance of two filters with grids crossed. 


Vertical Axis- wavenumber/wavelength 


Horizontal Axis- transmittance 
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APPENDIX B. 
SPECTRAL TRANSMITTANCE/REFLECTANCE OF THE 50 MM POLARIZER 


This appendix contains the Perkin-Elmer 
spectophotometer traces for the 50mm polarizing splitter. 
In Fig. Bl, the tracesmasseh the x designatecn are for the 
polarizer oriented 45 degrees to the incoming beam. The K,, 
trace is the combination of the traces and gives the overall 
transmittance of the splitter normal to the incoming beam. 
The xX symbols along the K,, trace are the overall 
transmittance of the splitter at 45 degrees to the incoming 
beam. It is seen that the peak occurs at a slightly shorter 
wavelength than at normal incidence. 

InJFig.B2, K,, Psmene vera li,ucm@ection efficiency in 


the configuration described in Chapt. VI. 
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APPENDIX C. 


MEASURED RADIANCE VALUES FOR THE NPS-NACIT SPLIT-FIELD 
POLARIME TER 


The following data were taken as input to the radiance 
analysis of polarizing splitter performance. 


fae Conditions: 


Room Temp 1/C, 

Source-Agema Black-Body @60C 

Polarized Source-BaFl IR Polarizer Pol%~99% 
Sensor-AGA 780 LW Band, f/1.87, Range 2 
Background Plate Temp- ~OC (heavy frost surface) 


eco ac 


fp 74.2 €, -11.2 
8.6 €, =5.4 
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APPENDIX D. 


CALIBRATION CURVES FOR FILTERED AND UNFILTERED AGA-780 
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Calibration for the AGA with no filters 


Figure D-1l. 
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Calibration Data for AGA with internal £ilteie 


Figure D2. 


al filters locatedurg 
Position 4 is the other 


A small difference is 


Pos 5 denotes one of the small intern 
position 5 of the AGA filter wheel. 
with orthogonal filter orientation. 


the two 


seen in thermal level with temperature between 


filters. 
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Figure D3. Calibratvemmeteve for tne Split-field 
Polarimeter. This shows level reading with temperature for 
the two channels of the split-field polarimeter. LWPOL is 
for the large BaFl external polarizer for comparison. Upper 
spot denotes the splitter transmitted channel imaged by the 
AGA. 
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APPENDIX E. 
COMPILATION OF MOSS LANDING IMAGES 
MAY, 1995 


The following is a complete image compilation of the entire 
Moss Landing experiment. The data tabulated is incomplete and 
only shows what is necessary for preliminary TBRD comparison of 
images. A sheet depicting a number of representitive images is 
appended to this list. Images corresponding to each of the table 
entries are available in digital form on hard disk at the NACIT 


MaD@ralOry. 


ane) 


The title header gives data about each of the images taken. 


Pol-polarizaeion Omientation of@the amage.gh,v-horiz,vert. x-image not 
(see image). 
Level-thermal level setting of the AGA control knob. 

Range-thermal level range setting of the AGA control. 
Offset-Offset of the calibration constant for image temperature compensation. 
Ovrflw, Undrflw-number of pixels in the image with radiance values above or 
below the 0-255 range. 
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APPENDIX F. IDL IMAGE PROCESSING CODES 


The following programs are used to determine TBRD 
automatically from a file which contains all the specific 
information about a target. The desired information is 
shown in the file compilation of image names in Appendix E. 

Also the IDL code to subtract the two split-images is 
contained in this appendix. The file headers denote the use 
of the code and the code itself is commented to provide its 


own detailed explanation. 


lez 5 


;fliral2.pro processes a single ship for contrast using sreadpol to get 
filenames from 

;experiment fliri2 

;Pete version that uses new contrast with no denominator 

;taken from flira8 which selects files directly 


72-17-95 revised x0,ix for x0 always to be peak. changed shipbox to find 
peak if the 
;new peak selected with the cursor is eliminated by the threshold 


;x0, yO peak value in contiguous area 

pale, LV, box coordinates 

;ship temperature ship 

; shipb raw data ship 

,yc horizon line point 

pAvZ Ly-Vve location cf Givisor in @ship Sox 


p RAE KKKKRKKEKEKKKKK KK KKK KKK KEK KEKE K EK KKK KKK KK KKKKEKKEKKKKEKKKKEKKEKKKKEKKKKKKKKK 
KHEERRKKKRKKKRKEKEKEKKRKEKK 


;routines readflir, shipbox, thresh, remake, fileonly,sreadpol used 
sreadflir reads file data from flir format 

FN @ ie Zon allows choice of horizon, formerly profiles 

7 shipbox finds contiguous ship and makes ship box 

sthreshld finds and returns edge threshold thresh from yy and t as well 
as peak x0, y0 

; remake uses threshold threshO to make ship and shipb from t and yy 
;fileonly separates filename from file and path for output 

;sreadpol reads polarization file and separates h or v in variable pol 


f!order=0 ;get rid of old program 
inverision 

iscale=2 ;display 2x size 
on_ioerror,problem ,if error, quit program 


PRK KKK KKEKEKEKEKEKREKEKEKKKKKK KKK KKK KK KKK KEK KEK KK KKK KKK KKEKEKEKK KEK KKKKKKEKEKKKKKKK KKK 


KREEKKEKKRKKKEKEKEKRKEKKEKKEKESE 


,select the files with wild cards 

pathi='' 

file=pathl 

if(n_elements (path) ne 0) then print, 'oldpath= ',path 

réedad, “Input Parwe ror £1 les or Creromeo|dpath',pathi 

if(pathi ne '') then path=pathit+'\' else path='.' 
if(strmid(path,strlen(path)-1,1) ne '\') then path=path+'\' 3 Ea 
oldpath error 

read, ‘input filename of polarization information', file 

sreadpol, file, files,pol , ;get files from flirts 
made separately 

files=path+files s;put files in standard 
findfile format 

filelist, files,path poecGram to lice 
filenames that are different 

7 LObe—O, Mune 1 G0 prince, ci les (1); sprint the filenames with 
a number 

y=! 

RESTART: read, 'Input start and stop number for processing CR to stop',i 
if strlen(i) eq O then goto, finish 

istart=f1x(1) 

istop=istart 

read, ‘Stop number = ",1Stop 


ZIG 


C5N=DBEARR ( (ISTOP-ISTART+1),3) ;The first z is full boxf™second, above 
horiz, third below 

read, 'input first guess difference t level, eg 1 to include all of 
sha. tdi f 

p RAR KKEKKKEKKKKEKEKKEKKKEKKKKKKKEKRKEKEKKKKK KKK KEK KKK KEK KEKKEKKEKRKKEKKKKKKKKKKKKKKKK 


KRKEKRKKEKKKEKRKEKEKEKEKEKK KE KK 


close,2 

fileout=' 

read, ‘input output filename for contrast results', fileout 
Seenw, 2, £1Lecut 

printf,2,systime(),' date of analysis' 

prmnt{t, 2mistep=dstart+1,' number of files' 


PRK RKKKKKKEEKKEKREKEK KKK KKKKEKKKE KKK KKKEKKKKEKKKKKKEKKKKEKEKKKKKKKKKKKKKKKKKKKKKK 


1 


KEKE KKK KEKE KEKEKEKEKEKEKEKEKEHE 


auto=' ; 

read, ‘input CR for semi-automatic proccessing', auto 

tdifsave=tdif ;save tdif so that changes will 
not affect batch 

printf,2,'Temperature difference = ',tdif ;save ina file 


meat, x,y, VC, 1¥ 
pmemeni,2,  Lilel, tdif,nx,ny, ye (horizon) ,iy(ship bottom) , Polarization’ 


Gee i-1istart, istop do begin ri = first of paar 

tdif=tdifsave ,restore the temperature to the 
Original batch value 

f=files (i) 

file=f ;save the name 

Beadtlinr,f,yy,c, dar, te, head ;read flir data, convert to 


temperature, add bar 
ee SR ARR R KK AKER ACT SDlay the raw 


Catark * KKK KKK KKEKKKKKEKKKKKKKKKEKKKKKKKKKKKKEKE 

loadct, 0 ;set color table O 

iscale=2 ;make pictures twice their normal 
size 

; remove the path length from the file string and label the tv 
image with the filename 

pathlen=strlen(path)+1 

filelen=strilen(file) 

labelen=filelen-pathlen 

label=strmid(file,pathlen, labelen) 
window, 0, xsize=140*iscale, ysize=140*iscale+50, xpos=0, ypos=0, title='0 
'+label 
tvscl, rebin(tt,140*iscale,165*iscale) ;display the image 
PERE KKEKKKKEKKEKKKKKKKKEKKEKKEEKKKKKKKEKKKKKKKKKKKKKE NG of temperature 
meee yA ee a x * KER x 


window, 2, xpos=0, ypos=165*iscale+35,xsize=140*iscale, ysize=140*iscale 


Meinct, ‘choose horizon with left button or reject picture with right’ 


tvscl, rebin(yy, 140*iscale,140*iscale) ;display the 
binary image raw 

HORIZON, smooth (rebin (yy, 140*iscale,140*iscale),5),xc,yc ;read cursor 
position 

if(!err eq 4) then goto, skiploop ;skip 


processing of file 
xc=xc/iscale 


yc=yc/iscale ;scale to 140 
size not picture size 
Print, "Xe,7C— mec, yo 


PRK RKKKEKKKEKKEKKEKKKKKKKKKKKKKKKECHreshold with temperature 
difference*** kK KK KEKKKKKKKKKKKKKKKKKKK 


ie] 


whreshld, twyy, thresO, x0, yO ,;define 
thres and find x0,y0 peak 


*RESHIP RESTORES THE VARIABLES MODIFED BY SHIPBOX, NOT THRESH 
errorct=0 sdeon’ t. allow intimate 


loopes of choices 
pRRRKRKKEKKEKEKEKEKKEKEKKEEKEKKEKKKEKKKKEKKRESHT DEKKKKEKKKKKKKEKKKRKKKKEKKKKKK 


KKEKEKKEKKEKKEKEKKKKEKEKKKKEKSE 


reship: 

remake, ship, shipb,t, yy, thres0O ;remake the 
thresholded ship and shipb using thres0O 

tvscl, rebin(shipb, 140*iscale,140*iscale) ;display the binary 
image raw 

tvers, x0*iscale, yO*iscale ,cursor on hot spot 
shipbox,shipb, ship, x0, yO, 1x, 1y,/ux,ny,tcceale 71X%,1ly,nx,ny dete 


ship box contiguous 
, WARNING ALL VARIABLES EXCEPT ISCALE ARE CHANGED IN THIS 


CALL 

boxdraw,ix*iscale,iy*iscale,nx*iscale,ny*iscale ;show where the box 
is in *iscale 

tvers, X0*iscale, yO*iscale, /device s;put the cursor on 


the peak location 
PRR KRKEKEKEKEKKRERERKEREKKEKKKEKKKEKKKKEREKEKKE EK KEKKEKKKKEKKEKEKKEEKEKKEKKEKEKEKEEKKKKEKK 


KKEKKEKKKKKKKKEKKKKEKKEKEKK 
;***above makes box with minimum 
thresholLGX** KKK KKKKKEKKEKKKEEKKKEKKKKKEKEKKEKEKEKKEKEKEKKKEEKEKEKEKEKE 
!'ERR=4 
zif(auto ne '* or errorct gt O) then begin suspend auto only 
for errors 
Peint, 'right button, use auto position, left, new Cursor position 
cursor, ixnew, iynew, /device, /down ;when pushed, 
ehoose 
print, '!ERR= ', !ERR 
if(!ERR eq 4) then print, 'auto selected' 
if(!'ERR eq 1) then begin 
xO=ixnew/iscale 
yO=iynew/iscale 


print, 'new peak at ',x0,y0 ;NEW peak 
at xO,yvo 
remake, ship, shipb,t, yy, thresO 7 LE new 


remade, the shipbox will have destroyed shipb 
*;WARNING SHIPBOX DESTROYS SHIPB AND SHIP 


shipbox, sht#pb, ship, x0, yO,1x,1y,nx,ny,iscale ;ixX,1lyY,nxXpe 
define ship box contiguous 
boxdraw,1ix*iscale,iy*iscale,nx*iscale,ny*iscale ;show where 
the box is in *iscale 
endif 
;endif ;semiauto is 


two button pushes now 


p RRRAAEKKEKXKKKKAKETE THE HOT SPOT IS FOUND DO LOCAL THRESHOLD IN THE 


BOX**XKRKKKKEKKKKKKK 


remake, ship, shipb,t,yy,thresOttdif ;remake with full 
threshold 

print, 'redisplay with full threshold’ 

Shipeosasnipbl, ship, <0, yO,1%,1lyY,;nx,ny,iscale 7x0, YO, nxpny 


define ship box contiguous 
tvscl, rebin(shipb, 140*iscale,140*iscale) 
boxdraw, ix*iscale,iy*iscale,nx*iscale,ny*iscale 


2S 


EVCES, 154 pplaee™ rie cursor on 
the max for a check 

print, 'cursor on) maximum 

pK EKEKEREKARE RE KER KAKAEEX*X TING a New Shipbox based venme prec e ecmerc 
CadLEXx RRR RRR KKKKKKEKKK 


if(total(shipb) eq 0) then begin LE data -isegall zero, 
threshold did not work 

errorct = errorerrs ;avoid an infinite 
loop 


print, 'no ship was found. Data is not uSabM@€ with this tdi£ft' 7 veit 
sbi (Ucbaielencere Kepe 7) ulgiaiel leyaxealiel 
Print,» Skil nee, ole 
goto, skiploop 
endif 
read, "input new tdsei*,; tait 
; auto=' ' suspend auto 
permanently on firsteerror 
goto, reship 


endif 

cai cf='a 

Pei, 'CaLrt= "7 tdit 

wat i— 

mealce ne '* or errorct gt 0) then $ 


mead, ‘input new tdif£ or CR for old or negative to skip’, tdiit 
if(strlen(tdiff) eq 0) then tdiff=tdif else begin 

tdif=float (tdiff) 

im(fcdi£ lt 0) then™gotopeskiploop 

goto, reship 
endelse 


° 
7 


PRK KEKKKKKEKKEKKKKKEKKEKKEKEKKEKKKEKEKKEEREEKEKKKEEKKKEKEKEKKEKKKEKKEKKEEKEEKEKKEKE 


KKREKKKKRKRKEKKKEKKEKRKEKEKEKE 


choices section checks to see which picture is larger and uses the 
appropriate tdif . 
;to remake and replot the temperatures reference windows 


;ship temperature ship in box & sv extension 
7; box temperature box not ship 

;shipb binary ship in box 

r;boxb binary box 


ship=(t gt thresO + tdif)*t ;be sure NOR to use 
the OLD threshold 
shipb=(t gt thresO + tdif)*yy 


Shoot (xs 1 xt, ys Ly ny) MSsubseer 
shipb=shipb (i1x:ixtnx,1iy:iytny) 
box=u eax tne ya yiny) ;everything in the 


box from original data 

boxb=yy (1x: ixtnx,iy:iytny) 

box=box* (ship eq 0) snot ship 
(temperature) in box 


wg 


boxb=boxb* (shipb eq 0) 
zscale=140/nx 
window, 5, xpos=140*iscale, ypos=0, xsize=nx* zscale, ysize=ny*zscale 


tvscl, congrid(ship,nx*zscale,ny*zscale) 


Z KREKKKRKEKEKKEKRKEKKKEKEKKKKKREKKKK KKK KE KKK KKK KK KKK KKK KKK KERR EKER KKK KKK KK KKK KK 


kaKw KKK KKKKK KKK KK KK KK 
~make CONtCrastsmmciay = first file, shipsudly cl Z=-abovewmhori zon, cile=up 


Gerhori zon 
,CSN full boxme above horizon, below horizon 


ve hOoELZonw ine Pole 
sly ship bottom 
7hz Vea iy, location of divisor in ship box 


. KKK KKEKKEKRKKKKKKEKKKE KK KRKEKKREKK KKK KR RKEK KK KEK KK KKK KKK KK KKK KEKE KKK KK KK KKK KK 
KRREKEKKKKEKKKRKEKKKKEEKEK 

hz=yc-1y 

pift(hz 1€C 0)) then horizon is belowsst@iipeand ship is flying 

,if(hz ge iytny) then horizon is completely above ship 


c5n(i-istart,0)=pcontra(shipb, boxb) 


PEEOMVeumnorlzony —-ly to ma 
if(hz le 0 or hz ge ny) then begin shz 0 means c5Sn (fz 
is ill defined as ship is on horizon 

Print, “hz enyeerwnz ny, (CrroL 7hz gt ny means 
c5n(1,1) is below horizon 

Gam (i-istartc, 1!) —con (1—-istart,0) 7hz 1t O means ship 


is flying above horizon 
c5n(i-istart, 2)=c5n(1i-istart,0) 

eneiinc sesp- oun 
Gomii=1stare, |) —pconerd (shilpb(~,hz:*)s boxb(* hz: *)) 
e5n (i-istare, 2) =pecntra(shipb(*, 0: hz voxb (+70: ha) 


;Guplicate the first if the value hz is zero 


endelse 

print, cplecnly (mies (1) raciw:, @—',tdif,nx,ny, yc, lv, * “~,pel (1) 
Print, CON ( (1-1SUaRT), 0), format="({(F7.3,2x), ‘Full ship’)” 
print, CoN((i-WsTART)),1),format—"( (F7.3,2x),'Above Horizon)" 
print, CoONse@-LSTARTieaZ), tormmata ((F7.3,2x),"Below Horizon*)™ 


printf, 2, fileonlemgGiskes( I eapaiehy), aa tdisipnx,ny, yc, iy," *",pol(® 
Drintt,2,CsnC(l-ISTART)., 0) Aroma ult). 3) 25), Fullesmip")" 

printf, 2, CSN((I-ISTART) ,1), format=" ( (EF 7 23hZx) , ‘Above Horizon 
printf, Z,€5N((I-ISTART) ,2), format=" ((F7.37@x) , “Below Howes zonvem 


skiploop: 


endfor 
ERR AEAKAKEEKEAKEEEREKEEERERKKERKKKKKEKKKEKEKKERERKEKKERERREREREE RRR KK KX He 


KKK KRKEKKKKKKKKEKKKEKKEKK 

Goto, Einish 

problem: print,'I had a problem and gquit' 
finish: print,dpmegram findshed' 

close,2 

end 


;program shipsz12 reads output of fliral2 (Pete Heisey) which is: 
;filename, tdif,ix,iy,nx,ny,H or V 


roo 


;and calculates TBRD with ship size and plots it 
close,1 shut up in case of an error 
tempp=' 


file=' ; 
read, 'input filename from analysis output',file 
Sn LOCEEOR, <COperr 
epenr,1,file 
nfiles=0 
1y=0 
Ve-1y 
NX=1y 
ny=1y 
i=0 
t=1.0 
Hv=' ! , EL Leteeune 
while not eof(1) do begin ;read to end of file 
which is not more than nfiles 
; loop: readf,1,tempp & print,tempp sread string and print 
nee 
loop: readf,1,tempp 
if(nfiles eq 0) then begin 7if number of files not 
moume yet, look for it 
if(strpos(tempp, 'number of files') gt 0) then begin 


reads,tempp,nfiles & print,tempp,nfiles get the number of files 

rdata=fltarr(5,nfiles) make array of number of 
files, may not use it all 

cdata=fltarr(nfiles, 3) yarray of contrast data 
; namearr=strarr(12,nfiles) array of filenames 


namearr=strarr(nfiles) 
pol=strarr(nfiles) 


print, 'finished making arrays, Nfiles= ' ,nfiles ;let me 
know arrays were made 
endif 
endif 
if(strpos(tempp,'.IMG') 1t 0) then goto, loop else $ 
namearr(i)=(strmid(tempp,0,strpos(tempp,'G')+1)) ;search for file end 
IMG 


sprint, (strmid(tempp, strpos (tempp, 'G')+1,strlen(tempp) )) 


y=strcompress(strmid(tempp, strpos (tempp, 'G')+1,strlen(tempp) ) ) 7get 
rid of white spaces 

meesGen ), Cyl X, Ve NXenVebV & DLINt, taux, 1V,nk,Ny pV ;convert data and 
Pent it 


MaEInttl, Zabeasip~nx, ngaayc,iye maee gthis 
Beads, V, tt; nx, away epaaa 
rdata(0,i)=t 
rdata(1,1)=nx 
rdata(2,1i)=ny 
rdata(3,1)=yc 
rdata(4,1)=1y 

pol (1)=HV 

y=1.0 

Beacdt,l,V & Cdata(i, 0) 
readf,l,y & cdata(i,1) 
readf,i,y & cdata(i,2) 
print, ‘incrementing loop 


ll 
axKNK 


ounter' 


eS 


1=itl ;increment the stored 


variable 
if(rdata(4,i-1) le 1) then i=i-1 s;if height is O, 
the data 

endwhile 

Goro, Scopp 

SUEODErr SPint, "1960 error’ 
goeo, loop 

stopp:close,1 

fe. tt 0) thengea—o0 

pDrEInt, 'fLinisheda,, 2 
rdata=rdata (*,0:1-1) 
cdata=cdata(0:1i-1,*) 
namearr=namearr(0:i-1) 
Prints Lije=—'y tile 

SPLOT Ye, RDATA, CDATA, FILE, pol 
end 


;program contrast.pro which calculates TBRD for each ship image 
;anayzed. it is a subroutine of shipsz12 

function contrast, ship, box 

contrast=(shipav (ship) -ship (box) ) 

return,contrast 

end 


;program splot 9 is subroutine of shipsz12 which takes the ship 
conatrast data and plots it vs. ship range determined by size 
PRO SPLOT9, RDATA, CDATA, file, pol 

ptitle=file 


Pedatat. 0) contrast for whole ship 

pecatat*,1) contrast for above horizon 

,ecdata(*,2Z) contrast for below horizon 

peaata(Z,*) ship size from height 

hpol=where(strtrim(pol,2) eq ‘h') sarray of number for h 
polarization 

vpol=where(strtrim(pol,2) eq 'v') ;array for v polarization 
help, hpol, vpol 


Lf (nmvelements (hpol) le 1) then hpol=where(strtrim(pol,2) née 'h") 


PsZ 


skip 


if(n_elements(vpol) le 1) themeeet—where(strerim (pol. 2 snow ae 
LOADCT, 39 

;read, 'Input title',ptitle 

;WINDOW, 0, XSIZE=1024, YSIZE=768 
;plot,1./rdata(2,*),cdata(*,2) , CHARSIZE=2, psym=2, coller=0, back@mouna—2Z55, 
charthick=2,title='!17'+ptitle $ 

;,Xtitle='Ship Range', ytitle='Target-Background Radiance 

Difference', /nodata 

POPlOt, lw/rdatal( 2) +) cdatede 2b pS VM=2.C01or—0 
Foplot,1./rdatea (Z,“), edatei 70) psym=1, coler=0 
poplot,1./rdata(2Z,*) ,cdata (*,1) ,psym=4,color=60 

;xyouts, .65, .87,charsize=3, charthick=2,'!17Above 

Horizon',color=60, /normal 
peaoucs,.65,.82,charsize=3,¢chareniek=2./417Full Shijel) cotor—0 -7normal 
;xyouts, .65, .77,charsize=3, charthick=2,'!17Below 

Horizon',color=60, /normal 

;WINDOW, 1,XSIZE=1024, YSIZE=768 
;plot,1./rdata(2,*),cdata(*,0O) , CHARSIZE=2,psym=2, color=0,background=255, 
eSWarthick=Z,title="!17'tptitle $ 

;,xtitle='Ship Range’, ytitle='Target-Background Radiance Difference’ 
poplom|./rdata(Z,*),cdata(*,1),psym=4,color=0 

ppoOplot, 1./rdata (Z,*) /cdata(*,2) ,psym=6, color=0 

WINDOW, 2, XSI ZE=1024, YSIZE=768 
plot,1./rdata(Z,*),cdata(*,2) , CHARSIZE=2, psym=2, color=0, background=255,c 
harthick=2,title='!17'+ptitle, /nodata $ 

,xtitle='Ship Range from Ship Size', ytitle='Target-Background Radiance 
Difference’ 

oplot,1./rdata(2,hpol),cdata(hpol,2),psym=2,color=0 
oplot,1./rdata(2,vpol),cdata(vpol,2),psym=6, color=0 
oplot,1./rdata(2,hpol),cdata (hpol,1),psym=1,color=0 
eelot,1./rdata(2,vpol),cdata(vpol,1),psym=4,color=0 

MyOuUCS, .65,.77,charsize=1, charthick=2,'!17 ** = Hee oy !17- 
V',color=0, /normal 

xyouts, .65, .87, charsize=1,charthick=2,'!17Above Horizon * 
'9B',color=0,/normal 

xyouts, .65, .82,charsize=1,charthick=2,'!17Below Horizon + 
19V',color=0,/normal 

end 


;this program subtracts the arrays of the two split images from the 
splitimage polarimeter. 
;program splitcom compares split image in yy array 


getfiles, files 
;get array of filenames in files 
EOD : 

;do it again 
read,'input number of file to get',ifile 
Mmotatile lt 0) then goto, problem 
f=files(ifile) 
readflir,f,yy,t,bar,tt,head Pde 
mre and Make temperature array 
!order=0 
On Crror, 1 
WeOt—yy (0. *, 703%) 

Veop-vyy(G: *7,0:69) 


Les 


yes= 7 

read,'input y to invert',yes 

if(yes eq 'y') then begin 

ybot=ybot (*, 69-indgen (70) ) ,invert the image 

ytop=ytop (*, 69-indgen (70) ) 

endif 

Hotspot, ytop, xhottmwhott 7FIND THE 

HOTSPOT IN THE IMAGE 

hotspot, ybot, xhotb, yhotb 7 SubrouTcrne 

window, 1, XSize=140*3, ysize=140*3 

GYVs¢el, Febiniveol, 1 40*37 70-3) 

tvscl,rebin(ytop, 140*3,70*3) ,0, 70-35 ;display both images 

xdif=xhott-xhotb 

calculate displacement to overlay hotspots 

ydaif=yhott-yhotb 

DEINE, KGa vai t.* sxdift -yait ' 

ytop=float(ytop ) make 

floating so subtraction will work 

ybot=float (ybot) 

displace,ybot, £ix (xdif) , faeetydl ; subroutine which 
;overlaps images 

window, 3, XSize=140*4, ysize=70*4 


absoldif=abs (ytop-ybot) ;make the 
values of the difference positive 

tvscl, rebin(absoldif,140*4,70*4) ;display difference 

SOCo,, Loop 


problem:print, 'exiting' 
;error and normal exit for -1 filename 
end 


;pro displace,yval,xdif,ydif 
ysize=size(yval) 


print, ysize 

if (ydif ge 0) then 

yval(0:*,abs (ydif) : ysize(2)-1)=yval(0:*,0:ysize(2)-1l-abs(ydif)) § 
else yval(0:*,O:ysize(2)-1-abs (ydif) )=yval (0:*,abs (ydif) : ysize(2)-1) 

end 


pro hotspot, yval,xhot,yhot 

ysize=intarr(5) 

ysize=size(yval) 

print,max(yval, ymax) & print, ymax 

print, ymax-(ymax/ysize(1))*ysize(1) ,ymax/ysize(1) 
xhot= ymax-(ymax/ysize(1))*ysize(1) 
yhot=ymax/ysize (1) 

end 
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APPENDIX G. 


SPLIT IMAGE POLARIMETER IMAGES, 3 AUGUST 1995 


The following appendix contains the images recorded of 
the USS John McCain (Arleigh Burke class DDG) on 3 August, 
1995 with the split-field polarimeter. This was a 
preliminary experiment conducted to test the performance of 
the optical system in the far field and to validate field 
performance of all system components. Any Data on the image 
file not contained in the Experimental notes is stored in 


the image file header. 


SD 


Experimental Notes 


Image File | Targe Thermal Remarks 
Name c Range 


faoeoaca [aa |__| setup test-amcblack boay source 
fmae02002 [aa | [infinite focus, ma @20" 
fmoeoaooa fae | i[? 
fmaeoa00e [ae | [black body teuperature s6deyrecs_ 
faaeoa00e |uccein [5 [| SSS 
jaoeoa007 |uccein [2 [SSS 
fmaeoa00e |uccain [2 [| SSS 
70803005 
mO803010 
mO803011 
70803022 


mois O0SOlS BB 


cal tab. 4004, BB 36 degrees 


m0803014 B 


B 
moseosoul> McCain 


moeos0i 9 FOV 





empty FOV with aft of ship on 









far right side 
mO0803020 McCain > 
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70803023 Focus set on inf, forbs 


iS i 
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